Scale-free networks without growth 
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In this letter, we proposed an ungrowing scale-free network model, wherein the total number 
of nodes is fixed and the evolution of network structure is driven by a rewiring process only. In 
spite of the idiographic form of G, by using a two-order master equation, we obtain the analytic 
solution of degree distribution in stable state of the network evolution under the condition that the 
selection probability G in rewiring process only depends on nodes' degrees. A particular kind of 
the present networks with G linearly correlated with degree is studied in detail. The analysis and 
simulations show that the degree distributions of these networks can varying from the Possion form 
to the power-law form with the decrease of a free parameter a, indicating the growth may not be a 
necessary condition of the self-organizaton of a network in a scale-free structure. 



PACS numbers: 



'.75.Hc, 64.60.Ak, 84.35.-K, 05.40.-a 



Many social, biological, and communication systems 
can be properly described as complex networks with 
nodes representing individuals and edges mimicking the 
interactions among them 0,013. Examples are numer- 
ous: these include the Internet, the World Wide Web, so- 
cial networks, metabolic networks, brain networks, food 
webs, and many others. Recent empirical studies indi- 
cate that the networks in various fields have some com- 
mon characteristics, the most important of which are 
called small- world effect Q and scale- free property 
The ubiquity of complex networks and the discoveries 
of the common network properties inspire scientists to 
construct general models. In the simplest way, these 
models can be categorized into two classes, the growing 
models and the ungrowing ones. In the growing mod- 
els, the number of nodes in network grows continuously, 
while in ungrowing cases, the number of nodes is fixed. 
The Watts-Strogatz (WS) model is the pioneer and rep- 
resentative one of ungrowing models, which can be con- 
structed by starting with a regular network and randomly 
moving one endpoint of each edge with probability p |4| • 
WS networks display small-world effect, that is, they are 
highly clustered and of small average distance. The most 
successful one of growing models is the Barabasi- Albert 
(BA) model @, which suggests that two main ingredi- 
ents of the self-organization of a network in a scale-free 
structure are growth and preferential attachment. These 
points to the facts that most networks grow continuously 
by adding new nodes, which are preferentially attached 
to existing nodes with a large number of neighbors. 

By using the mean-field theory on a toy model, 
Barabasi et al assert that the growth is one of the nec- 
essaryconditions for the emergence of scale-free prop- 
erty [6|, and the networks without growth may display 
an exponential degree distribution. This hypothesis is 
now widely accepted. However, some real-life networks 
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having fixed size (or almost fixed size) are of scale-free 
property. For instance, consider the friendship networks 
of school children wherein the students are represented 
by nodes and two nodes are connected by an edge if the 
corresponding students have good friendship Q>[g. These 
networks do not grow with time since the number of stu- 
dents in a class is almost fixed. Some students have many 
good friends, while many others have only a few friends. 
The students' heterogeneous societal ability leads to a 
skewed degree distribution (or a skewed strength distri- 
bution for weighted case 8]), approximated to a power- 
law form. A typical example in biological systems is 
the functional networks of human brain @. Although 
its structure varies with time, the brain functional net- 
work's sizes is unchanged [To| while it displays scale-free 
property. Another significant examples are the relation- 
ship networks of countries like world trade web 0, 
and world exchange arrangements web 0, These 
networks are scale-free, while their sizes almost do not 
change in the recent a few years. 

Although these networks mentioned above have (al- 
most) fixed size, they are not static since their struc- 
tures vary with time. These motions can be consid- 
ered as rewiring processes, that is, some existing edges 
are removed while some new ones are added. For ex- 
ample, in friendship networks, good friends may quarrel 
about beliefs, money, or some other things, and become 
impassive to each other; while some ones will become 
good friends for their common interests and difficulties. 
The economic globalization promotes different countries' 
economies to be integrated together in terms of trade 
and capital flow, making countries more interdependent 
than ever. In order to avoid financial crises, more flex- 
ible economic regime should be implemented, and the 
government may change the exchange-rate/trade regime 
for certain economic and political reasons, leading redis- 
tribution of degree/ weight in world exchange arrange- 
ment networks and world trade networks. These redis- 
tribution/rewiring processes make the structures vary- 
ing ceaselessly in recent years, but the scale-free prop- 
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Let P(k, t) represents the number of nodes with de- 
gree k at the time step t. P(k, t) satisfies the following 
normalization conditions 

^P(k,t) = N, (1) 



Y,kP(k,t) = 2E, (2) 



J2G(k)P(k,t) = l. (3) 

fc 

It is straightforward to write down the master equation 
for P(k,t), 



FIG. 1: (Color online) The degree distribution of the present 
model for the case a 3> 1. The black hollow squares and 
red solid circles represent the analytic and simulation results, 
respectively. The analytic result agrees with the simulation 
accurately, and both of them obey Possion forms. 

erty is always observed. Notice that, in the ungrowing 
model proposed by Barabasi et al @, no edges will be 
removed thus no rewiring processes will occur. In this 
letter, for the first time to our knowledge, a scale-free 
network model without growth is proposed, and we argue 
that the ungrowing networks can have power-law degree 
distributions, which attributes to the rewiring processes. 

For simplicity, we consider a network model with both 
fixed number of nodes N and fixed number of edges E. 
The initial state of the network can be a random graph 
or any other types of graph. Then the network evolves 
based on the following rewiring processes: At each time 
step, an edge is randomly selected and removed from the 
network. At the same time, a node is selected with the 
preferential probability G(fci), where ki is the degree of 
the ith node. Another node is selected also with the 
above preferential probability. Then a new edge con- 
necting these two nodes is created. Notice that in this 
rewiring process, the total number of edges is unchanged. 
The above process is repeated at each time step. Finally, 
we expect that the network reaches an equilibrium state 
which is independent of the initial state. 



, . (k + l)P(k + 1, t) - kP(k, t) 
P(k, t + 1) - P(k, t) = y - >-A — _lJ Lj_2 

+2G(k - l)P(k -l,t)- 2G(k)P(k, t). 

(4) 

It may be helpful to explain the physical meaning of var- 
ious terms in the above equation. The first two terms 
on the right hand side represent the net gain of P(k, t) 
due to the edge removing process. In particular, the first 
term represents the net gain of P(k, t) when the remov- 
ing edge connects a node with the degree k + 1. Notice 
that when an edge is randomly selected, the connecting 
nodes have more chance to be of large degree. Explicitly, 
the node with degree k is selected with the probability 
kj J^k kP(k, t) — k/2E. Since each edge connects to two 
nodes, the first term is thus (k + l)P(k + l,i)/E. Sim- 
ilarly, the second term represents the net loss of P{k, t) 
when the removing edge connects a node with the degree 
k. The third and fourth terms represents the net gain of 
P(k, t) due to the edge adding process. The third term 
represents the net gain of P(k,t) when the adding edge 
connects a node with the degree k — 1. Notice that a 
node with degree k is selected with the probability G(k). 
The fourth term represents the net loss of P(k,t) when 
the adding edge connects a node with the degree k. 

When t is sufficiently large, we expect that P(k, t) ap- 
proaches a stationary state denoted by P s (k) that satis- 
fies 



p, + m(h + d - kP , (k) + G{k _ m{k _ t) _ g(m{k) _ ^ (5) 



I 

where P s (k) also satisfy the normalization conditions order equation in which P s {k + 1) is determined by P s (k) 
Eqs.(l-3). Notice that this master equation is a second and P s (k — 1). Eqs. (1-3,5) can be solved by the following 
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method. Define 



and 



H(k) = G(k)P a (k) 



W{k) = kP s {k)/2E. 
Then Eqs. (2-3,5) can be rewritten as 

£V(*) = i, 

k 

J2H(k) = l, 

k 

H(k) = H(k - 1) + W(k + 1) - W(k). 
By iteration, one can obtain 

H(k) = H(0) + W(k + 1) - W(l). 
Since W(0) = and 

£fr(fc) = XV(fc) = l, 

k k 

one immediately have 

W(k + 1) = ff(fc) 



(6) 
(7) 
(8) 
(9) 
(10) 
(11) 



(12) 



for k > 0. Substituting Eqs. (6-7) into Eq.(12), one ob- 
tains the following solution 

P.(k) = ^p G ( fc - l) G ( fc - 2)...G(0)P S (0) (13) 

for k > 1 with P s (0) determined by the normalization 
condition Eq.(l). It should be stressed that given the 
dependence of G(k) on fc, there is still a proportional 
coefficient in G(k) determined by Eq.(3). 

Hereinafter, we focus on a special form, also one of the 
simplest cases, of G(k), where G{k) is linearly correlated 
with k, 



G{k) 



k 



2E + Na 



(14) 



where fej denotes the degree of node i, and a is a con- 
stant (Notice that G(k) satisfies the normalization con- 
dition Eq. (3).). Substituting the above equation into 
Eq. (13), one can straightforwardly obtain, in this spe- 
cial case, that 



P s (0) = N{- 



a + 7 



-r 



(15) 



Ps(k) = N( 
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FIG. 2: The degree distribution of the present model with 
a — 0.01 and 7 = 5. The hollow squares and solid curve 
represent the simulation and analytical results, respectively. 
The analytical result agrees with the simulation well, obeying 
an approximately power-law form. The inset exhibits the case 
of a = 0.01 and 7 = 10 for comparison. 



where 7 = 2E/N. 

When a ^> 1, the mechanism of preferential attach- 
ment is destroyed, that is, since a probably plays the 
main part in the sum a + ki, the probability a high degree 
node being selected is approximately the same as that of a 
low degree node (see Eq. 14). Therefore, one expects that 
the present network will approach a completely random 
network with its degree distribution obeying the Possion 
distribution [131 



Ps(k)=N— 7 fc . 



(17) 



When a <C 1 and 7 « 1, most nodes are isolated and 
the nodes with large degree are very few. However, when 
a <C 1 and 7^1, interesting result emerges. The cor- 
responding degree distribution obeys a power-law form 
with exponent approximated to 1 when k <C 7/a 



P.(k) = N(-)°j. 

7 k 



(18) 



(16) 



In succession, we will show some simulations. 10 6 time 
steps arc performed in all simulations while only the final 
2 x 10 5 time steps are used for statistical average. The 
initial graph is simply the completely random graph with 
given N and E. The network size if fixed as TV = 1000 
unless a special statement is addressed. In figure 1, we 
report the case for very large a, where p s (k) := P s (k)/N 
denotes the normalized probability function for degree 
distribution. The simulation result strongly supports the 
analytic one (see Eq. (16)), and both of them obey Pos- 
sion forms. In this figure, one will find that there are 
about 5% nodes with degree zero. However, this proba- 
bility only reflects the average number of isolated nodes 
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FIG. 3: (Color online) The degree distributions of the present 
model. The parameter 7 = 10 is fixed, and the black squares 
and red circles represent the cases with a = 1.0 and a — 
0.1, respectively. The degree distribution with a — 1.0 obeys 
approximately an exponential form rather than a power-law 
form. 

over all configurations. Because of the rewiring mech- 
anism, no nodes are always isolated. Although real-life 
networks are not always connected, most previous stud- 
ies have focused on connected graphs. Therefore, here- 
inafter, we neglect the isolated nodes. 

In figure 2, we report the analytical and simulation re- 
sults about degree distributions generated by very small 
a = 0.01 and 7 = 5. The degree distribution displays 
obviously scale-free property with a cut-off at k » 120. 
The solid curve denotes the analytical solution, which 
agrees with the simulation well. We also have checked 
that the departure from Eq.(18) will get greater if the 
parameter a becomes larger. In the inset, the case of 
7 = 10 is shown for comparison. As mentioned above, 
the two extreme cases with a w and a 3> 1 will lead 



to the power-law and Possion distributions, respectively. 
In figure 3, we report the simulation results for a = 0.1 
and a = 1.0. The departure from the power law becomes 
larger as the increase of a, and the degree distribution 
with a = 1.0 obeys approximately an exponential form 
rather than a power-law form. 

In conclusion, we proposed an ungrowing model and 
obtained the analytic results for arbitrary forms of selec- 
tion probability G(k). One special kind of the present 
networks is investigated in detail, which can generate 
networks with different classes of degree distributions by 
tuning two parameters. Especially, when a is close to 
zero and 7 M, the present networks display scale-free 
property with its exponent approximately equal to that 
of world exchange arrangement web [l4| . 

The previous empirical studies have demonstrated the 
extensive existence of scale-free networks. Most of them 
are growing at all times, while some others are of (al- 
most) fixed size. Most previous models Q 0, 0, IE @ 
suggested the growth a key ingredients of the emergence 
of scale-free structure, thus fail to provide an underly- 
ing mechanism by which the ungrowing networks exhibit- 
ing scale-free property. Here we argue that the rewiring 
mechanism widely existed in real-life networks like friend- 
ship networks [71 la , brain functional networks |9|, world 
trade web [Iull2| . and world exchange arrangements web 
[l3l FLU , may be the origin of the emergence of scale-free 
structure in fixed-size networks. 

Finally, we would like to point out that many different 
classes of degree distributions can be obtained from dif- 
ferent forms of G(k) |16j . In particular, one can obtain 
the power law degree distribution P s (k) ~ /c~ /3 ~ 1 with 

tunable exponent (3 by setting G(k) ~ [Eip wnen the 
number of edges E is properly chosen |16| . 
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